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ABSTRACT 


With the growth in demand of electricity for vari- 
ous purposes, today power system has become complex in 
nature with a number of interconnections of different ge- 
nerating stations located at different places and supplying 
composite loads connected to it at different voltage le- 
vels, The generating stations include both hydro as well 
as thermal plants connected in a grid. In such a system, 
if the installed capacity exceeds total demand plus losses, 
a power engineer will have a number of strategies to meet 
the load demand, but the best will be the one which results 
in the minimum operating cost satisfying all the operating 
constraints on the system. Several programming techniques 
have been used in the past to get an* optimal solution. The 
present work as reported in this thesis exploits the tech- 
nique of dynamic programming to arrive at a realistic answer 
in solving a short range scheduling problem of a hydro-ther- 
mal system including losses. Achievements of the present 
work are increased accuracy, faster rate of convergence, 
reduction in computational efforts and simple algorithm with 
an easy mathematical approach. 
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CHAPTER I 


INTRODUCTION 


In the present day when economy at every walk of life 
is a basic requirement, the electric consumers will always look 
forward to the cheapest electrical energy made available to 
them* For a power system to return a reasonable__ profit on 
the capital invested, various modes of operation exist right 
from generation to distribution to arrive at an economic one; i 
but of all, economic load scheduling is important and challengj 
ing* Rates fixed by regulatory bodies and the importance of 
conservation of fuel place extreme pressure on power companies 
to achieve maximum efficiency of operation and to improve it 
continuously in order to. maintain a reasonable relation betweei 
cost of a Kwh to a consumer and the cost to the company of 
delivering a Kwh in the face of constantly rising prices of 
fuel, labour, supplies and maintenance. In operating the sys- 
tem for any load condition, the contribution from each plant 
must be determined so that the cost of delivered power is a mil 
mum. This is the challenge before an engineer and the central 
theme of the present work is to meet this challenge limiting 
the considerations to the hydro-thermal system on a short^rang 


basis* 


2 


In the earlier days, there used to be mostly small 
thermal plants for each locality known as urban systems with the 
sources close to the load centre* But with the considerable . 
growth- in the demand of electricity and rising cost of fuel and 
its limited reserve, a need was felt to include hydel plants 
also alongwith the thermal plants in a grid system with th% 
basic aim to achieve minimum operating cost and exploit the 
capability of hydel plants having higher reliability and greater 
speed of response when run as a peak load plants However, to 
achieve these, a number of constraints have to be satisfied. In 
an interconnected system, the generators may not be close to the 
load centre and hence, transmission losses will occur which have 
to be met to arrive at a realistic optimal solution resulting in 
minimum operating cost. The operating cost of thermal plants 
and also that of hydro-thermal plants is mainly the cost of fuel 
as the cost of water is negligible* Cost of labour, maintenance 
etc. are assumed to vary as a fixed percentage of fuel cost. 
Hence, cost function of hydro-thermal plant is defined as a 
nonlinear function of plant generations. For a short range prob- 
lem, the level of water in the reservoir of hydel plant remains 
almost constant for a given time interval. Only constraint 
will be, therefore, the specified quantity of water to be uti- 
lised in each hydel plant in the given time interval. By eco- 
nomic or optimal load scheduling we mean to determine the gene- 
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rations of different plants, hydel as well as thermal, such 
that the total operating cost is optimal and various constraints 
as mentioned above are also satisfied. The problem, therefore, 
is the non-linear programming problem having nonlinear equality- 
constraints. 

Break through in the area of optimal generations of hy- 
dro-thermal system was probably first achieved by C3iandler, 
Dandeno, Glimn and Kirchmayer They used Lagranges multi- 

plier to convert the nonlinear objective function with the non- 
linear equality constraints into an unconstrained objective 
function* The optimum scheduling problem as considered here is 
defined as follows: 


While using specified amount of water from the hydro 
plants over a given period of time, it is desired to minimize 
the. total fuel cost (in dollars). If it is assumed that, for 
the period of time under consideration, the hydro plants operate 
at essentially constant head, then the cordination equations 
whose solution gives the optimum schedule are 


dF, 


n 


dL. 


dP 


4* X 


sn 


sn 


= X ; n=l , 2 , 


a 


( 1 * 1 ) 


d w • ^ L 

d^fc dPHj 




= X» 


j=0H-l,a+2, ... 


( 1 . 2 ) 
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Where, 


dFn 

-Twr* = incremental fuel cost of thermal plant n in dollars 
sn 

per Mwh 


^ W . ’ 

= incremental transmission loss of thermal plant n 
sn 

dw. 

dPH ' ~ incremental water discharge rate at hydro plant j 
3 

. in cubic feet per Mwh 


dL^ 

’ 6 ' PH " ~ incremental transmission loss of hydro plant j 


\ = incremental cost of received power in dollars per 

Mwh 


= constants which convert an incremental water dis- 
charge rate into an equivalent incremental plant 
cost. 


a - number of thermal plants ; w. = water input to 
plant j in cubic feet per second 


p-a = number of hydro 
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For the purpose of evaluating the savings possible by 
scheduling generation according to solution of co-ordination 
equations, the following methods of scheduling had been con- 
sidered: 


!• The Co-ordination Equations Method : 


i.e. if the (ciF„)/(clP-_) and (dw.)/(dPH.) terms in 
equations (1*1) and (1.2) are assumed linear functions of out- 
put power, they can be written in slope intercept form as 


dF 

dPTT “^nn^sn‘*’^n 


(1.3) 


and 

dw. 


(1.4) 


dF dw. 

By substituting the values of (^^^ ) and (’ ^ ' ^" ) from 

s n j 

equation (1*3) and (1.4) into equations (1.1 ) and (1.2) res- 
pectively, the choice of \ then determines a set of generations 

P and PH. and thus a particular received load, 
sn j 


2, The Equal Incremental Plant Costs Method : 

If the terms (dL^)/(dPg^) and (dL^/dPH^) in equations 
(1.1 ) and (1.2) are neglected, then again substituting the 



6 


dF dw . 

value of “tp — and 735^ from (1,3) and (1*4) into (l.l) and 
sn 

(1.2), we get, 

‘’nn^'sti'^^n + yWj = \ 

Again, the choice of \ determines a set of generations 
Pg and PH. and thus a particular load. 

3, The maximum Efficiency Hydro Method : 

It is a proven practice to operate hydro plants near the 
■point of maximum efficiency except during times of surplus 
water. 


While working with these coordination equations on digi- 

[21 

tal computers, Dandeno observed that the direct application 
of these co-ordination equations give solutions which sometimes 
dictates generations outside the plant capacity and also nega- 
tive generations of certain plants since restrictions regarding 
the plant capacities had not been included in the problem for- 
mulation. 

However, significant improvement was achieved when 
Wije Perera^^^ and others applied Pontryagin* s Maximum Prin- 
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ciple followed by Kuhn-Tucker conditions to solve the problem 
imposing also the upper and lower bounds on the operating 
range of each plant in the form of inequality constraints. 

Nonlinear programming technique was used by Sokkappa 

[ 4 ] 

B.G. • Here transmission losses were included by the use 
of loss coefficients. For each subinterval of optimisation, 
the constraint that is most likely to be violated, was picked 
up and a slack variable was associated with the constraint* 

The gradient of the objective function i.e. cost of the gene- 
ration , was evaluated and the steepest descent method was 
used to obtain the solution of the problem, starting from a 
known initial solution i.e. generation schedule. However, 
this method was found to be inefficient when the system size 
becomes large. 

fsl 

Gopala Rao*- solved a short range problem using Lasdon*s 
decomposition technique^ splitting the problem of higher di- 
mensions into subproblem of smaller dimensions to be solved 
iteratively. The method suffers from a drawback i.e. the ini- 
tial choice of the dual variable vector, which plays a signifi- 
cant role in the overall computational time, is crucial. A 
brief account of the method is given below. 
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The optimisation interval is subdivided into *K' equal 
intervals such that the load demand during the interval at all 
the stations remains constant. To find transmission losses » 
A.C. load flow equations were used. The variables in the prob- 
lem are bus voltage magnitudes and their angles. Maximum and 
minimum limits on active and reactive power generations and 
maximum limits on bus voltage magnitudes are imposed. 


Let a system has 'M' hydro and 'N* thermal stations and 
*m’ total buses. Then, during interval, the injection 
equations can be expressed as. 




#) = 


m . 
E 

j=l 




pK Y 


cos 


'Sij -«1 




( 1 . 5 ) 


Ji(6^EK) = 


- ^ ^ 


( 1 . 6 ). 


where , 

K th 

E = bus voltage magnitude in K interval 

K th 

6 = bus voltage angle in K interval 

Y. ./9. .= ij^^ element of Y-bus matrix 

+w 

I., J. = injections into the system at i bus. 


Therefore, power generation equations can be written as 
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Qf(6^E‘^) = Ji(6^E'^) + 

for i =. 1, .••; (M+N) 
and K = 1 , « • • > K 

whore and are active and reactive power generations at 
bus, and and are active and reactive load demands 

•f’H 

at the same bus in the K interval. 


Similarly, the load demand equations will be. 


I, ( + Poi = 0 j 


1^(6'^, E*^) + Qci = 0 


for i = M+N+1 ••• m 


and K=l, ...,K 


The objective function, i.e. the cost of generation 
can be expressed as. 


f(6'',E'') = I ''2'' C,(P^(6^E^)) 

K=1 i=M+l 


(1.7) 


Now, the optimisation problem can be mathematically 


stated as 



Minimise f(6^,E^) 

subject to the following constraints! 

(i) Load Demand Equations, 

+ PdI = ° 

Jj(6'^,e’^) + = 0 

for i = M+N+1, m 

and K = 1, , . , , K 

(ii) operating limits on the system. 


pinin < 




< 


Q?(6^E^) <0?®* 


pK , pmax 

ti V 


for i = 1, •••» M+N 
and K = 1, •«*, K 


(iii) restrictions to use specified quantity of water, 
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where, is the volume of water available for power generation 
The objective function can, therefore, be written as, 


f(6^,E^) = Z fj, 

K=1 ^ 

where, f^(6^,E^) = Z C.(P^(6 ,E^)) 

i=M+l ^ ^ 


(1.9) 


Associating dual variables Uj^,**«, with constraint 
equation (1»8), the Lagrangian function can bo written as, 


L(6^,E^,u) 


Z fj,(6^,E^) + 
K=1 ^ 


M K 

Z u. ( Z 
i=l ^ K=1 


yi(pf(6^,E^)) - ) 

( 1 . 10 ) 


Therefore, we have a dual problem, i.e. 

Maximise h(u) over all u 

K K \ 

where h(u) = rain L (6 ,E ,u) 

^ for K = K 

The set also, known as constraining set to which the 
sub problem variables 6^ ^ ^9 be constrained is 
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defined by constraints (j.) ^ (ii) mentioned above* 

Using Lasdon’s decomposition technique^^^ , the Lagran- 
gian function is minimised by solying the following 'K* sub 
problems, for K = 1, , K. 

i.e* minimise fj^(6^,E^) + 2 Uj^y^(P^(6^,E^) ? 

subject to the constraints (i) and (ii) mentioned above. 

Therefore, we can say that the crux of the 
method mentioned above is to have an additively separable 
nonlinear programming problem split into subproblems and solved 
iteratively* The selection of dual variable is extremely 
crucial as the minimisation process is accomplished by it. 

Ill 

Hano, Taraura Y and Harita^S ^ have used continuous 

fisl 

maximum principle of L.S. Pontryagin •' to the most economic 
operation of Hydro- thermal systems. A brief account of the 
method employed by them is given below. 

Let us consider, the system having two plants j one 
hydel and the other thermal, jointly supplying electric power 
to a load centre through a system having lossless transmission 
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lines. The rate of water in flow into the reservoir as J is 
assumed to be constant. Load demand is represented by the 
daily load curve spanning over 24 hours. The fuel cost of 
thermal plant is approximiated by the square function of the 
thermal output [i.e F(G) = a^^G + a^G ]. 


Let the water storage in the reservoir, the rate of 
water discharge, and rate of water inflow be expressed as x^Ct), 
u(t) and J respectively, then 



It is assumed that the effect of evaporation,, leakage, 
or overflow due to rainfall etc. have not been taken into 
account in equation (1.11 ) above. Turbine lower and upper bounds 
rates of water discharge are given by 


u < u < u 


( 1 . 12 ) 


The power output P of the hydro plant may be expressed 
approximately as 


p = Hq (1 + C Xj^) (u - q) 


(1.13) 
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Where, = basic water head 

o 

C = correction factor of water head 
q =s rate of non-effective water discharge 

The power balance equation may be expressed as, 

= G + P (1.14) 

where, Pj^ = the load demand 

P = Power output of hydel plant 
G = Output of thermal plant 

Here, transmission losses have been neglected. 

With these definitions, the problem of the economical 
load dispatching may be initiated in the simplest form as 
follows: 

Given a time-dependent function of load as Pp^(t), and 
the initial and final storages of the water as and Xj^(T), 

respectively, determine a control u such that the integrated 
fuel cost of the thermal plant 

T 

4) = / F(G)dt (1.15) 

0 
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be minimal, where F(G) is the fuel cost* 

Employing the maximum principle theorem, the problem 
will be restated as follows: Given a time dependent function 
Pf^Ct), the initial and final states of the system as x(0) « 
[xj^CO),©] and x(T) = [xj^(T),T], respectively and a set of diffe- 
rential equations governing the dynamics of the system as. 




o 

dt 

= F(G) 

(1.16) 

dXj^ 

= J - u 

(1.17) 

dt 

dX2 

— 1 

(1.18) 

dt 

— -L, 


determine the optimal control u which minimizes the integral 
function. 

T 

x^(T) j F(G)dt (1.19) 

The optimizing condition of the maximum principle states, 
that, the optimal control here, is the one which maximizes the 
corresponding Hamiltonian function 

H = - F(G) + (J - u) + ^2 


( 1 . 20 ) 
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where, and 1 IJ 2 are the adjoint variables to the system and 
are governed by a set of differential equations i.e. 


dlll^ 


= - i 

(1.21) 

dt 


dlllj 

dt 

_ -dH . 

dXj ■ 

dF 

“ “ dG dt 

(1.22) 


and. also, ” ^1 “ dG ^o^^ ^^1^ ” ^1 

(1.23) 


making (1.23) equal to zero, we get 
dF _ 

dG ■“ Hq( 1 + Cxj_) 


(1.24) 


Since, G is a function of u as is evident by equation 
(1.13) above, knowing from equation (1.24), G can be ex- 
pressed as. 


G = 



‘•H.U+xJ 



(1.25) 


If G lies outside its maximum and minimum values G and G, 
modifications in the choice of u is made to bring back G with- 
in the permissible region. From equation (l.i7) and (1.19) we 
get. 


I 
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u < u <. u 

u > U 
u < U . 

Since the maximum principle approach involves numerical 
integration of differential equations * 

i = 1,2,, ..,n ; ^ ; i = *•*» 

terval for each step of numerical integration should be small 
enough depanding upon the degree of accuracy desired. As per 
the authors view, dynamic programming was found to be less effi 
cient in terms of computation time and accuracy when compared to 
their proposed method. However, an effective method has not yet 
been developed to determine the initial values of the adjoint 
variables to a system given by equation (1.20) above. This is 
the serious drawback of the system and it has been accepted by 

the authors also, 

Srivastava K.N. and Singh have used the successive 

dual linear programming technique to obtain solutions to the 
complex problem of optimal scheduling of hydro-thermal systems 
including transmission losses for a short range problem. It has 
been claimed that the technique has succeeded in reduction of 



R 



(1+Cx^) 

“• T* q XI 

- u 


if 

= u 


if 
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storage and computational efforts. A brief description of the 
mathematical formulation of the work is presented below* 

Let the cost characteristics of different thermal plants 
and discharge characteristics of hydel plants as a function of 
plant generations are represented by, 


and 




YiCh^) = + Qj^h^ + x^(h^) 


(1.26) 

(1.27) 


where, 



i 


3 

K 


y,(h^) 


= Generation of thermal plant j in MW 
= Generation of hydel plant i in MW 
- 1,*.., Mi.e. M hydel plants 

= 1,..., Ni.e.N thermal plants 

= Sub intervals of optimization (usually 
a day descretized into K equal inter- 
vals ) 

= cost function of j thermal plant in 
Rs/hr 

= the discharge function of i^^ hydel 
plant in m /hr 




c.»Pj^,q 


i'^i 


are constants* 
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While considering the cost associated with a particular 

K 
3 


interval K, the subscript K is dropped from S. and h^ and the 


objective function is written for that particular interval as 
follows : 

v N M n 

L (a,+b.s.+c.sp + E Y. (P 4 +q.h.+r.hf) (1.28) 

J u J u J i=l ^ ^ 11 11 


subject to constraints 
N M 


(a) 

E S. + E h. = P^ + P, (S.,h. ) 

j=l J i=l ^ ^ 3 1 

(1.29) 

(b) 

eHiin X c / cHiax 

(1.30) 

and (c) 

hf" < h, < 

1—1—1 

(1.31) 


where i 


Pq = demand in MW on the system 
= transmission losses 


Linearizing the scaler cost function of equation (1.28) 
by Taylor series expansion around the initial operating state 
(Sj,h°), the incremental cost function is obtained as, 

n r. ^ o 

= I (b. + 2C. S?) <fcS. + E Y4(q4+2r.h?) a h. 
j=l J 3 3 3 1 1 1 1 1 


AfK = 
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1«Q. 


Af^ = 


M+N 

Z d. aP, 
i=l ^ ^ 


(1.32) 


where = incremental cost co-efficient 

AP^ s= incremental change in the generation 
of plants. 

The constraints of the equations (1.30), (1,31) and 
(1.29) can be rewritten as 


g. = S. - > 0 


(1.33) 


^N+i = H > 0 for i = 1, ..., N 


(1.34) 


_ v,max 
’2N+M+i “i 


h^ 0 for i = 1, • • • , 


(1.35) 


IN m 

32N+2M+1 “ ^i “ ^D,“ ^ 


(1.36) 


Linearizing the equations (1.33) to (1.36) around the 
current operating state by Taylor series expansion, the new 
linearized constraints are obtained as follows. 


= aS^ > - S? for i = 1, ..., N (l.37a) 
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%+i = - > S° - sf ^ for i = 

^2N+i ~ ^ ~ ^ ~ 

^2N+M = “ t for 1 


1 ,..., N (1.37b) 
1,..., M (1,38a) 
= 1,...,M (1.38b) 


92N+2MH-1 = V^L^S^h?) - I S° 


M 

- Z 
i=l 


(1.39^,) 


Equation (1.32), together with equation (l.37a) to (1.39_) 
forms a linear procrarraning problem in incremental variables 
A S . and Ah.. 


Using Dantzig’s method to make problem variables to be 

non-negative, a non-negative vector Z of dimension (N+M) with 

elements is defined as follows, 

-aS^ - ASj 

for j = 1, . , . ,N 

(1.40a) 

^N+i ” A h^ -Ah^ 

for i = 1, . . . ,M 

(1.40b) 

where, A =’ - 

sP 

(1,41a) 

A hf ^ = hf^ . 

v,° 

’ ^i 

(1.41b) 



22 


substitution of equations (1,41a) and (l,41b) in the 
expression for incremental objective function given by equation 
(1,32) yields, 

ginin 

Af^ = d’^Z + d"^ [ ^min ] (1,42) 

In equation (1,42), since second term is constant, the 
objective function reduces to, 

Af^ = d^Z (1.43) 

Then the LP problem can be stated as : Determine Z 
that minimizes the objective function of equation (1,43) subject 
to the constraints of equations (1,37) to (1.39) after substitut- 
ing the values from equation (1.41a) and (1.41b), 

A comparison has been made with the present work with the 
above mentioned dual LP technique taking a common problem, 
which shows that a greater accuracy has been achieved in terms 
of cost per unit supplied using dynamic programming technique 
as used in the present work (as reported in this thesis), 

(91 

Singh L.P, and Aggarwal used Dynamic programming 

technique to solve a short range problem where transmission 
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losses were not included to arrive at the optimal solution. 
A brief mathematical formulation of the above work is given 
below. 


The problem is to determine the scheduling of genera- 
tion of the hydro-thermal plants i.e. to determine S. and 

sJ 

h^ such that 

N- M 

Cj ^ Z C.(s.) + I Yiy4(h.) (1.44) 

^ j.=l ^ ^ i=l ^ ^ ^ 

subject to the constraints 

N M 


(a) 

21 S- 4* 21 h. = Pr 
j=i j i=i ^ 

+ PT(S,h) 

(1.45) 

(b) 

5”?“ < S. < S?®* 1 
3 - J •- 3 3 

1 

(1.46) 

and hf < h, < h"'?^ 1 

^ 4. JL ^ 



(c) 

/ Yiyi(h.)dt = 


(1.47) 


Assuming that a total number of R plants which satisfy 
the total demand D, the plant supplies the demand 2* 
Therefore, R-1 plants will supply demand D-Z with the minimum 
cost of generation Fp^_2^(D-Z). Hence using, principle of opti- 
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mality, we .get, 


Fp^(D) Min [Uj^(Z) + Fp^^j^(D-Z) 

where, = cost of generating Z by the plant, 

in general, - 

[Uj^(Z) + Fj^__^(D-Z)] where, N=number of 

thcfrmal plants. 

To include hydel plant to run as a peak load plant, we 

put 

R = N+1 i.e. first hydel plant is put into operation when 
all N thermal plants must be operating; the cost function we 
get, 




Min 


h: 


mm 


1 Z 5 a 


max 


[Tiyi(hi)+FN(D-Z)] 


*i ~ -^“i 

of course, checking that equation (l.4l) is satisfied. 


OUTLINE OF THE THESIS 

The objective of the present work is to obtain solutions 
to the complex problem of optimal scheduling in hydro-thermal 
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power systems including transmission losses using the techni- 
que of Dynamic programming. Merits of this scheme are increased 
accuracy, faster rate of convergence, and simple algorithm. A 
brief outline of the work reported in the different chapters 
is given below. 

The second chapter introduces the mathematical theory 
in brief of the technique of dynamic programming due to 
R. Bellman^^^^ and others' 

The third chapter describes a general approach to eco- 
nomic load scheduling of Hydro-thermal systems due to Kirch- 
mayer^^"^^, Singh Stevenson Jr. 

In the chapter 4 the formulation and solution procedure 
of the problem of optimal load scheduling of hydro-thermal 
systems using the technique of Dynamic programming have been 
described. A case study comprising of three thermal and one 
hydel plants has been presented to illustrate the solution 
procedure. Transmission losses have been considered through 
the use of loss coefficients. A comparison with successive 
dual LP technique in solving a common problem has also been 
presented which affirms that increased accuracy, faster rate 
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of convergence, and simple algorithm are the assets of the 
present work (reported in this thesis). 

Finally, chapter 5 concludes with the overall summary 
of the work reported in this thesis. A few proposals for 
further scope of investigations in this area are also made at 
the end. 



CHAPTER II 


DYNAMIC PROGRAMMING 


2.1 INTRODUCTION 

*Dynaniic Programming’ is a mathematical tool for sol- 
ving problems in the catagory referred as multistage decision 
processes^^^^ , The name was so given because Dynamic Progra- 
mming in many cases leads to a solution in the form of a 
program for a digital computer. Dynamic Programming was de- 
veloped during 1950 's by R. Bellman^ to study the optimi- 
zation problem arising in the industry, economics and in the 
social services where optimization technique of linear and 
non linear programming and calculus of variation and its ge- 
neralization are not applicable. Such catagory of problems 
has been referred to as multistage decision process. Econo- 
mic load scheduling of hydro-thermal system essentially belongs 
to the catagory of problems known as the multistage decision 
process requiring a sequence of decisions and which can be 
tackled more effectively by the technique of dynamic progra- 
mming, 

2.2 MULTISTAGE DECISION PROCESS 

To understand the multistage decision process, let us 
consider a system which can be described at discrete times 
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by a finite number of variables (known as state variables) 
as The values of these state variables at time 

(say stage) i is denoted by X 2 (i)» x^(i)». In the 

form of state vectors we can write, 

X = [x^.Xj,. 

and its value at stage i is given by 

X(i) = [xj^d), X 2 (i), ... x^(i)] 

Now, at stage i = 1,2,..,N, we can make a decision 
d(i) from a number of possible choices. The effect of taking 
decision d(i) at stage i will change the state of the system 
i.e. at stage i+1, the state vector is X(i+1). 

At stage i, of the system a return function R will be 
associated which is of the form, 

R(i) =R Cx(i), d(i), i] 

i.e. the return function R depends upon the state of the 
system at the stage i.e. X;(i) decision d(i), and stage i. 
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Restrictions can also be imposed on the return function 
to take integral values only and which can be of the form, 

0 [X(i), d(i), i] < = € 

Because, there are N stages and at each stage if a 
single decision is made then there will be sequence of decisions 
such as d{l), d(2),.,,d(N) and corresponding to each such 
decisions , there will be a sequence of returns. The total 
return I, for the whole process may be expressed as, 

I = R(i) + R(2) + R(N) 

N 

or simply I = Z R(i) (2.1) 

i=l 

The above system is called Discrete Time Multistage 
Decision Process. 

Now we will take an example of multistage decision 
process. The example concerns an investor, who, each year 
must allocate a sum of money between two possibilities} either 
reinvest it or spend it. The problem is specified as follows: 

(i) Let Rs X bo the amount of money available in any 
year, 

(ii) Let Rs U be the amount of money allocated to in- 
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vestment, so that Rs(X-U) is spent in the same 
year, 

(iii) Let us assume that allocation decision is made only 
once a year and that years are numbered 1,2,.., i. 
Then the amount of money available in year i is 
written as Rs X(i), the amount of money invested 
in year i is Rs U(i) and the amount spent in 
year i is Rs [x(i) -- u(i)], 

(iv) Assume a constant interest K (K>1) so that as a 
result of investing Rs u(i) in year i the amount 
of money available, for allocation in the next year, 
i.e. i+1 is x(i+l) = Ku(i), 

(v) Let H(i) be the amount of satisfaction derived from 
spending money in any year i, 

(vi) Assume that the satisfaction H(i) depends only on 
[X(i) - u(i)] i.e. the amount spent in year i, 
and that the relation between satisfaction and 
amount spent does not change with time and can be 
described by a known function of [X(i) - u(i)]as 

H(i) = H[x(i) - u(i)] 

This is called the return function, 

(vii) Therefore, we can say that the total return I, 

derived by the investor from N consecutive years 
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i = is given by, 

I = H[x(l) - u(l)] + H[x(2') - u(2)]+ •..hCx(N) - u(N)j 

( 2 . 2 ) 

, N 

or simply i = Z H(i) (2,3) 

i=l 

Thus I depends on the initial amount of money x(l) 
and on the series of decisions u(l),u(2),. .u(N), 

(viii) Each decision u(i) is subject to restrictions because 
the amount of money invested in any year cannot be 
greater than the amount of money available that year 
and cannot be less than zero, i,e. 

0 < u(i) < x(i) j i =:1,..,N • 

This is an example of multistage decision process. 

2,3 PRINCIPLE OF OPTIMALITY 

The principle of optimality, orginally as stated by 
BELLMAN^^^^ is applicable to all multistage decision process 
and is to the effect that : 

*An optimal policy has the property that, whatever be 
the initial state and the initial be the first decision d(l), 
the remaining decisions constitute an optimal policy with regard 
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to the state resulting from the first decision’. 


Therefore, our problem is to find the optimal value of 
I from equation (2,1), 


i,e, if the initial state is S and d(l)>d(2). ,,d(N) 
decisions are taken at stages i = 1,2,,,,,N, Then, optimal 
return over N stages vail be expressed as. 




max 

d(l),d(2),,,d(N) 



(2,4) 


2.4 DYNAMIC PROGRAMMING TECHNIQUE 

The technique of dynamic programming depends upon the 
principle of optimality which is applicable to a system wrtiere 
the state of the system at the stage i+1 depends only upon 
the state of the system at the previous stage i, decision d(i) 
taken at stage i and the stage i itself. 

The general form of the transformation in a determinis- 
tic multistage decision process satisfying the above condition 
is represented by N deterministic functions and can be 


written as. 
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* ^^[xj^{i)i ••• Xj^(i ) jd( 1 ) j il 


It is convenient to use vector notation by defining a 
vector function G having components Gj^ ••• i.e* 


G2 

j 

G « 


% 


and to represent the transformation by 


x(i+l) « G[x(i),d(i),i] 


and the total return function can be written as 

N ■ 

I == 2 R(i) 

i=l 

The problem is to find a sequence of N decisions, one 
decision at each stage for which the value of I is maximum* 
The same technique is precisely used in the formulation and 
solution of the problem envisaged in the present work which 
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will be described in Chapter iv later. Let the initial state 
of the system is 5 such that, 

fl^(s) = Maximum value of I when N stages remain 

Let us take the first decision d(l) which is arbitrary^’, 
then corresponding to this decision the state of the system 
becomes » 

x(2) = G[s,d(l)] 

and there remains N-1 stages to go. The principle of 
optimality states, that, when the optimal policy is used, the 
return from the remaining N-l stages starting from state, 
G[s,d(l)] is 

fN-l[G {s,d(l)} ] 

But the return from the first stage is, 

R[ 6 ,d(l )3 

So the total return is, 

I = R[s,d(l)3 + (s^dd^j (2.5) 

When the optimal policy is used, this total return I is 
maximized with respect to all N decisions, including the first 
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decision d(l), so the optimal return over the N stage is 

f^is) « ^ (s,d(l)} ]) (2*6) 

subject to specified restrictions 

or ~ R[s,,dj^(s)3 + [G[s,dj^(s)) ] (2#7) 

Where dj^(s) is the first optimal decision to be taken 
when the system starts with the state S and total N stages to 
go. 

Equation (2.7) is a functional recurrence equation from 
which the optimal return functions can be found. The 

initial solution can be found for a single-stage process, by 

t 

taking N=l, when the definitions give 

max 

f,(s)* [R(s,d(l))] 

^ d(l) ^ 

= R[s,d(l)] 


Similarly, from equation (2.7), putting N=2, gives, 


f^Cs) = R[s,d2(s)3 + f^ [G (.s,d2(s)} ] 


( 2 , 8 ) 
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Equation (2,8) can be used to compute f2 (s)» continuing 
like this, we can compute fjq(s) and dj^(s) for any value of N . 
Here, gives us the maximum value of I for the N-stage 

process when the initial state is s. 

From the above discussion, it is evident that we have 
to solve a maximization (or minimization) problem in one vari- 
able only at every iteration. This is infact, the great force 
achieved by the use of dynamic programming. This technique 
has been used to solve the problem of optimal load scheduling 
of hydro-thermal system including transmission losses which 
has been brought out in chapter 4, 



CHAPTER III 


ECONOMIC LOAD SCHEDULING OF HyPROTHERMAL SYSTEMS 


3.1 INTRODUCTION 

The Economic Load Scheduling of Hydro-Thermal Systems is 
a complex problem and it can be solved by considering it basicatlly 
a long range or a short range problem. The long range problem 
normally takes into account the constraints spanning over a 
year. The inflows into the reservoirs show an annual cycli- 
city. There may be seasonal variations in power demand on the 
system. Raihfall and head variation of reservoirs from season 
to season also play an important role in detemining hydro- 
generation. Various methods have appeared in the literatures 
for determining economic load scheduling of Hydro-thermal 
systems taking long range problems into consideration. Since, 
the present work is a short range problem, it is beyond the 
scope of this chapter to explain long range problems. 

The work reported here is a short range problem having 
an optimisation interval of day consisting of 24 hours and 
load demand on the system has been considered for subintervals 
of one hour each. The solution to the problem will stipulate 
the quantity of water to be utilised over a day. This quan- 
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tity of water must lie within the specified quantity of water 
over the period (a day)* Therefore, the solution to the short 
range problem consists of an optimal plan for utilisation of 
this v/ater for power generation and the corresponding optimal 
set of thermal generations, considering a load demand on the 
system over a day (24 hours), of course satisfying various 
constraints imposed on the operation of the system. 

Several attempts have been made in the past to solve a 
short range economic load scheduling problem of hydro-thermal 
systems. A brief overview of different techniques used in the 
past has already been given in chapter 1 with a brief intro- 
duction to their mathematical approach. However, a general 
approach using Lagrangian multipliers as described by Kirch- 
mayer^^^^, Singh L.P*^^^^, Stevenson^^^^ is presented here in 
order to get acquainted with the requirements, constraints and 
mathematical explanation of a short range hydro-thermal system* 

3,2 LAGRANGIAN. MULTIPLIER APPROACH FOR SOLUTION OF ECONOMIC 
LOAD SCHEDULING OF HYDRO-THERMAL SYSTEM 

The main cost of operation of hydro-thermal system is 
the cost of fuel for the thermal plants, because it has been 
assumed that the problem is of short range type and hence. 
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there will not be any appreciable change in the level of water 
in the reservoir during the interval under consideration , 
However, specified quantiy of water must be utilised with- 
in the interval at each hydel plant i. 


Thus, based upon above assumptions, the problem is to 
find the thermal generation and the hydel generation h^; 
both of which are functions of time such that the total cost 
function (i.e. objective function) 


T N , , 

C-r = / E C.(S.)dt (3.1) 

i 0 j=l J ^ 

is minimum subject to the constraints, 

(a) E S. + £ h. = Pt + PT(S,h) (3.2) 

j=l J i=l ^ ' 


(b) 


where. 


c.(s.) 


S 


0 ^ ^ 


(3.3) 


fuel cost in Rs/hr of the j formal plant 

3 "tin 

Turbine discharge in m per hour of the i 

hydel plant 

S S«,...»S»^ = Generation of N thermal plants 
12 H 
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h - ^2^ • # • = Generation of M hydol plants 

= Demand in MW of at any time t 
Pj = Transmission losses as a function of plant gene- 
ration 

T = Final time defining the short range problem 
under study* 

The equation (3.1) is a non linear objective function 
with both the constraint equations (3.2) and (3.3) being non- 
linear, such problems are referred as constrained nonlinear 
programming problem. The nonlinear objective function with 
the nonlinear equality constraints can be convered into an 
unconstrained objective function F by proper choice of Lag- 
ranges multiplier X and y 

T 

i.e. F = - y^if y^(hj^)dt - Y^) - 

N M 

( E S. + E h. - P. - Pt) (3.4) 

j=l ^ i=l ^ ^ ^ 

Using law of calculus for minimization, the above funct- 
ion F can be solved by taking partial derivatives of F with 
respect to the plant generation S and h and equating them to 
zero for minimum i.e. 



41 


6F 


6S. 


dC. 

° “ dS^ 


6Pj 


4E- = n ^ 




0 = Y 


dYi 


i dh- 


+ \ 


6P^ 

6tr 


-X 


or. 


dC. 

"dgr 


6S. 


ili. 

dh^ 


6Pn 


+ X 


6h^ 


= X 


(3.5) 

(3.6) 


(3.7) 


(3.8) 


The direct application of the above co-ordination equa- 
tions (equations 3.7 and 3.8) results in solutions which 
sometimes dictates plant generations outside the plant capa- 
city and also negative generations of certain plants. This 

f2l 

phenomenon was first observed by Dandeno^ while working with 
the co-ordination equations on the digital computer . 


Therefore, the problem is extended by imposing upper and 
lower bounds on the operating range of each plant as indica- 
ted by the inequalities (equations 3.9 and 3.10) 


< S. < for j = 1, N 

3 - J - J 


h?^'' < h. < h?®^ for i = 1 

1 -■ X “* 1 


M 


(3.9) 

(3.10) 


f « • • > 
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Where and are the minimum and and 

ni3X 

^i maximum limits of operation of the corresponding 

plant. Dandeno also doubted about the constancy of y^ if 
these limits on the operating range of plants are included. 

[3l 

Wije Perera made an attempt to solve the problem by app- 
lying Pontiyagin*s maximum principle, the problem thus ob- 
tained is to maximize the Hamiltonian H at each instant 'of 
time i.e. 

N M 

Maximize H- - £ C.(S.) - E Y 4 y^(h.) (3.11) 

j==l J J i=l ^ ^ ^ 

The value of should be so chosen that the condition 
given by the equation (3.3) is satisfied because known as 

3 

water value, actually converts the water consumption in m per 
hour to the cost in Rs per hour. An increase in the value of 
results in lesser water consumption at the i hydel plant and 
vice versa* Our aim is to minimise the function given by the 
equation (3.11). 

Since the minimization of any function is equivalent 
to the negative of the maximization of the same function, 
therefore, minimization of the cost function Cj is given by 
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Minimise C ^ = - H 


N M . , 

= Z C. (S.) + Z YiYiCh^) (3.12) 

3=1 ^ 3 i=l 


The minimisation of the cost function Cj as defined by 
equation (3.12) corresponds to the minimization of the total 
operating cost of thermal as well as hydel plants. This is 
because the first terra in equation (3.12) represents operating 
cost of the thermal plant and the second term, that of hydel 
plant. 



OIAPTER IV 


OPTIMAL LO AD SCHEDULING OF HYDRO-THERMAL SYSTEM 

INCLUDING LOSSES 

4.1 INTRODUCTION 

The ultimate aim of load scheduling problem is to mini- 
mise the cost of generation of power while satisfying all the 
constraints imposed on the system. Several techniques inclu- 
ding Non-linear programming, Pontryagin's maximum principle, 
successive dual linear programming and dynamic programming etc. 
were used in the past to solve a problem of short-range hydro- 
thermal systems [1,2,3, 4, 5, 7, 8, 9]. Each has its merits and 
demerits* 

In the present work, a short range problem of load sche- 
duling of hydro-thermal system including transmission losses 
has been taken for study and solution. This chapter is devo- 
ted to the problem formulation using dynamic programming 
technique, its solution procedure and development of simple 
algorithm and finally testing it for a given test system. In 
a combined hydro-thermal system, almost negligible cost of 
generation is associated with hydro stations as the cost of 
water is negligible* However, there is always a constraint 
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in the form of availability of water (neglecting the probability 
of failure of plants) i,e» constraint on total energy to be 
generated over a day. 

4.2 PROBLEM FORMULATION USING DYI'^IC PROGRAMMING 

The optimal load scheduling of hydro-thermal system 
primarily belongs to the catagories of problems known as multi 
stage decision process and, therefore, is quite suitable for 
formulation and study using the technique of dynamic progra- 
mming^ In the case of optimal load scheduling of the 

hydro- thermal system, we are required to take a sequence of 
decisions and at each decision point a number of strategies 
exists. Ihe decision that we are required to take in this case, 
is regarding the generation schedule of different plants. Here, 
the stage corresponds to the number of plants that are in ope- 
ration and the decision at any stage corresponds to the proper 
choice of generation schedule. Changing the decision at any 
stage will bring about a change in the state of the system. 
Therefore, the effect of taking decisions for generation sche- 
dule for a particular plant will automatically bring about a 
change in demand to be met by the remaining plants. With this 
change in the demand i.e. change in the state, everytime a 
return function in the terms of cost function of the hydro- 
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thermal system will be obtained} out of which one of the cost 
functions will surely be optimal because of taking a sequence 
of decisions regarding the generation schedule. 

The total number of units both thermal as well as hydel 
available, their individual cost characteristics and the load 
cycle on the station are assumed to be known in advance. It 
is also assumed, that, the load on each unit or combination 
of units changes in small but uniform step size. We can 
start arbitrarily with any two units. The most economical 
combination is determined for all the discrete load levels of 
the combined output of the two units. At a particular load 
level, the most economical solution may be to run either of 
the units or both the units with a certain load sharing bet“* 
ween the two units. The third unit is then added and the above 
procedure is repeated to find the most economic answer. The 
process is repeated till all available units are put into 
operation. The advantage of dynamic programming (DP) is that, 
having obtained optimal way of loading M units, it is quite 
easy to determine the optimal manner of loading N+1 units. 


The aim of applying DP is to know the optimal scheduling 
of hydro-thermal plants. That is to determine and h^ such 
that, total cost function Cj given by. 
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N 

^ C.(S.) + 

j=l J J 


M 

Z 

1=1 


Yiyi(hi) 


(4.1) 


is optimal subject to the following constraints 


(a) 

N 

M 


E 

j=l 

S. + E h. = P,+P^(S,h) 

^ i=l ^ ^ ^ 

(4.2) 

(b) 

(.min 

< Sj < Sj 

(4.3) 

(c) 

-min 

*'1 

< < hf =< 

(4.4) 

(d) 

T 



/ yi(hi)dt = 

(4.5) 


It is obvious that the optimal value of Cj depends 
upon the total number of stages i.e. total number of plants 
to be considered (M+N in this case) and the initial state of 
the system i.e. initial demand D which takes into account the 
constraint equations (4.2), (4.3) and (4.4). 


Thus ^ 

Let Fj,(D) = Minimum value of (equation 4.1) when 
there are K stages that is K plants in operation and the 
initial state is D (i.e. initial demand is D). 
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Let us assume that out of a total number of K plants 
which meet total demand D, plant meets the demand G* 

This is the first decision we have taken. The effect of this 
decision will be to meet demand D—G by the remaining K— 1 
plants. Let the cost of generating G by the plant be 
fj^(G) such that, 

- Thj,(G) } if plant is thermal 

and fj^(G) = Y^YkCG); if plant is hydel 

Therefore, 

^K-l (D-G) = Minimum cost of generating D-G demand from 
K-1 plants which are in operation, of course having a check 
that demand D-G lies between the sura of the lower bounds and 
the sum of the upper bounds of K-l plants. Using the principle 
of optimality* we get the following relation, 

Min 

Fk(D) = [fj^(G) + Fj^^j^(D-G)3 (4.6) 

G 

Here, it is to be noted that consequent to the initial 

ju w 

decision that K plant supplies arbitrary demand G, Fj^__,j^(D-G) 
corresponds to the remaining K-l decisions which are always 
optimal with regards to the state resulting from the initial 
decisions i.e. first decision which is arbitrary. 
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According to formulation, when the first thermal plant 
is put into operation, we get, by putting K=l, 

min 

Fj^Cd) = f (D) = f (D) = C,(D) (4.6) 

G 


^min 

^1 


Here, demand D must satisfy the constraint equation 
< D < 


To operate the second plant also, we put = 2 and 
from equation (4.6), we get, 

, ^ Min Min 

F2(D) = [f2(G) + Fj^(D-G)] = [C^iG) + C3^(D-G)] 

^ (4.7) 


Similarly, by putting K=3,4.... etc, we go on finding 
F 2 (D), F 4 (D) ...... etc till all the plants are put into ope- 
ration. 


While initiating hydel plant, a proper value (i.e. 
water value) of is assigned and at the same time, the hydel 
generation also is maintained in the limits i.e. h£^'^<,h^jCh?^^ * 
The value of the total discharge obtained in the case of hydel 
plant is compared with the specified quantity of water to be 
utilised by this hydel plant in the specified interval i.e. 
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T 

/ Yiyi(h^)dt = Yj (4.8) 

Initially, a low value of y is assigned and its value 
is incremented till equation (4*8) is satisfied and also 
demand D is met 

N M 

i.e. D < ( Z S. + Z h. ) 

j=l ^ i=l ^ 

Thus, a generation schedule is obtained without including 
transmission losses. It is not possible to assess these losses 
in the beginning as the plant generation (loading) upon which 
these losses depend are not known in the beginning. Therefore, 
problem remains to obtain an optimal generation schedule inclu— 
ding transmission losses. The same has been done in the follow- 
ing paragraphs. A solution procedure is also suggested. 

4.3 SOLUTION PROCEDURE TO INaUDE TRANSMISSION LOSSES USIN G 
dynamic PROGRAMMI^j(g 

Use of transmission losses as function of plant genera- 
tions is made to derive the following transmission loss formu- 

ri4l 

la as given by Kirchmayer*’ Here, in the loss formula, 

the linear term and the constant term have been excluded. And 
hence, transmission loss is given by. 


Equation 


where, P 


and P 


'‘tec. ^ 


88321 


*44 1011^0 


51 



N 

L 

n=:l 


P B ^ 
m mn n 


(4.9) 


where, P^ = Transmission loss 

Pjjj,Pn = real power generation at m, 
n^^ plants 

= loss coefficients also known 
mn 

as B coefficients which are 
constant under certain assumed 
operating conditions such as 
unchanged location of genera- 
tors, particular transmission 
network and particular plants* 

.9) may be written in the matrix form as, 

Pl = P^BP 






h' 


®11 

®12 • • * • 

®1N 

^2 


®21 

^22 * ♦ « • 

®2N 

# 


# 



• 

and B = 

• 



♦ 

1 

# 





®N1 


CQ 





J 


is transpose of P» 
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Using dynamic programming, a generation schedule is 
obtained without including transmission losses in the first 
stage of optimisation* Let this schedule be P • In the second 
stage, transmission losses using B coefficients are calculated 
which are given by, 

(4.10) 


Now, lot us assume that P changes slightly to include 
transmission losses for an optimal solution. From equation 
(4.10) above, the incremental transmission loss (ITL) at plant 
n may be expressed by. 


dP 


L 

n 


2ZB P 
mn m 


(4.11) 


Let the cost of generation is Fj^ for a particular load 
level D. Since, transmission losses are generally less than 
15-205^ of total load supplied, we can assume that individual 
plant generation changes only slightly while optimising to 
include losses* As a result, the per unit cost of power F^/D 
and the partial derivatives reasonably be assumed 

to be constant in the second stage of optimisation for a 
particular load level. 

Per unit cost of power = Fq/D 
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F dP 

Change in the cost of transmission w.r.t. P =(^)(3p^) 

n u 

(4,12) 

Let us assume that P +a P be the optimum solution re- 
quired to include transmission losses. Assuming quadratic 
plant cost curves and witting it for the n^^ plant as shown, 
we got 


^n ^n ^In^n ^2nE" 


(4.13) 


If P^ changes to (P^ +a P^^), then. 


(4.14) 


or, 


4C„ = (C^^ + 2P^Co«)a P„ + a P 


n 


n"2n'“ "n ’ '"2n “ "n 


(4.15) 


Now, for including cost of transmission, we will modify 
incremental cost function as follows, 


dP, 


^ ^n “ ^^In ^^n^n^^ ^n ^2n'^ ^n (4.16) 


Incremental Load =aL = P, 


(4.17) 


For, the incremental lower and upper bounds, the following steps 
are followed, 

(i) If Pjj = 0 ; , A UB^ * 0 
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where, a UB^ and are the incremental upper and 

lower bounds of plant n 

(ii) If P^?^0 

we consider UB^ - P and 0.15 P (i.e. 15?^ P_) 

n n n n 

A UB^ = minimum of UB_-P_ and 0,15 P 
n n n n 

similarly, a LB^ = maximum of (-0.15 P^, LB^ - P^^) 

(iii) Having found a LB^^ and ^UB^ for all plants in the 
schedule, the following check is applied 

if Z aUB < Pt , then step (ii) is repeated with 0,2 P^^ 
n " 

and so on till the following condition is satisfied 

Ea UB > Pr 
n L 

Having determined a a UB^, incremental cost 

function from equation (4.16) and incremental load from equa- 
tion (4.17), use of dynamic programming is again made to find 


The final generation schedule is given by ^ 


for a load demand D. 
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4.4 CASE STUDY 

A short range load scheduling problem of 24 hours duration 
has been considered and the period has been subdivided into 
interval of one hour each. The demand at the generating stat- 
ion is assumed to be constant during each hourly interval and 
at the end of each interval, the demand increases or decreases 
in jumps or remains constant. The system considered consists 
of throe thermal plants and one hydel plant. The cost chara- 
cteristics of thermal plants and discharge characteristics of 
hydel plant "are given in Table I* 

The transmission losses have been calculated for the 
given system by taking the assumed values of loss coefficients 
as shown in Table II. 

The above problem has been solved on DEC-10 Digital Com- 
puter at IIT, KANPUR using the formulation and the solution 
procedure discussed in the proceeding sections for the load 
demand shown in Table III. 

The optimal generation schedule giving optimal cost of 
generation, individual plant generations and transmission losses 
for a particular load level is shown in Table III. Solution to 
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the problem has been obtained for a water value ( X ) of 0.131* 

A flow chart as shown in Fig. 1 has been attached in the appendix 
based upon which a computer programme was developed, 

A comparison has been made with the solutions obtained 
by the technique of successive dual linear programming. It 
has been found that greater accuracy has been achieved by 
using the technique of dynamic programming* 


TABLE - I 

Cos t/dis charge Characteristics 


Lower Upper 

Bound bound cost characteristics 
in MW in MW 


(a) Thermal plant 


50 

,200 

C^=100+0.1 S^+0»01 ] 


- 

60 

170 

C 2 = 120 + 0.1 S 2 + 0»02 

Rs/hr 

- 

50 

215 

0^=150+0.2 S 3 +O.OI S 3 i 


- 



(b) Hydel plant 



15 

65 

* 2 

y ^=140+20h j+0 . 06h 


500.00 


cubic feet/hour 


Quantity 
of water 
to be uti- 
lized over 
24 hours 
in 
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, TABLE - II 



Loss 

Coefficients 


0.0005 

0.00005 

,0.0002 

0,00003 

0.00005 

0.00004 

0.00018 

-0.00011 

0.0002 

0.00018 

0.0005 

-0.00012 

0.00003 

-0.00011 

-0.00012 

0,00023 


\ 
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CHAPTER V 


CONCLUSION 

In this chapter a broad review of the work reported in 
thLs fch-esis and a few proposals for the scope of future in— 
v^sidgia tLofis in this area have been presented* 

Utie optimal load scheduling of hydro-thermal system 
is a ccmplex problem. Several attempts were made in the past 
to sol"ve jsuch problem including nonlinear programming techni- 
ques, ffliaxEiouKtt principle of Pontryagin and successive dual 
ILncair purogrramming. Each has its own merits and demerits. 

Since*, e-con oniic load scheduling of hydro-thermal systems belongs 
to the category of problem known as multistage decision process, 
the t eohJiique of dynamic programming which is an appropriate 
tna thiGinaitdcall tool for such categories of problems has been 
usiet dm -thio present work for a realistic optimal solution. 

In this to iMiulation, transmission losses were also included by 
the use of loss-coefficients. Load demands spanning over a 
period of 24 hours have been taken into consideration and 
those hare been treated as deterministic. A case study has 
bean pres ented by taking a short range scheduling problem 
horiny th lee thermal and one hydel plant. The algorithm de- 
veloped is very simple and it does not require any complex 
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mathGinatical approach* This is infact an edge over various 
other techniques used in the past, to tackle a short-range 
problem. Choice of water value y is very important to solve 
tho problems of this nature* However, by incorporating 
Binary search in the present work to find ^ considerable 
reduction in solution time is achieved even with flat start 
of say 1,0 and 0.0 respectively in the beginning. 

The solution to the present problem was found for a water 
value of 0.131. If we start from y = 0.0 in the beginning 
and adopt sequential search, we have to go through as many as 
131 iterations to come to the solution of the problem. But 
by using Binary search we have to carry out only log 2 131 i.e. 
about 8 iterations to arrive at the final solution which is 
certainly a great improvement over sequential search. The 
final generation schedule is more realistic and optimal cost 
is found to be more accurate and acceptable. A comparison 
has been made with the solutions obtained by the technique 
of successive dual linear programming. The proposed method 
is very effective fot a system which is predominantly thermal 
with a few hydel stations. In case the number of hydel plants 
is more, they can be replaced by an equivalent hydel plant 
which will have the total energy equal, to the sum of the 
energy of the respective hydel plants provided these plants 
are located at about the same distance from the load centre. 
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SCOPE OF FUTURE im/ESTIGATIONS IN THE AREA 

Further investigations in the area are proposed conse- 
quent to the achievements made in the present work. They are 
as follows: 

1. In the above problem, a deterministic load demand has 
been assumed but in practice, there may be uncertaint- 
ies in the load demand i.e. making the problem a proba- 
bilistic one. Also discharge of water may be uncertain 
over the period under study, it is, therefore, proposed 
to make suitable modifications of the solution procedure 
to take into account these uncertainties. 

2. In the formulation of the problem, if reactive power 
also is taken into consideration, further savings can 
be achieved. Because, both real and reactive power . 
contribute to the transmission losses. In that case, 
the linear and the constant term of the loss formula 
also are to be taken into account. Use of load flow 
analysis could be made to determine reactive powers and 
the line flows and they can be included in the problem 
as additional constraints. But at the same time, it is 
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quite obvious that computational time will increase 
and therefore, it would be important to chalk out a 
balance between the computational efforts and time 
required and the optimal answers obtained. Although, 
according to Kirchmayer^^^.^ it was found that the 
additional savings obtained by the use of loss-formula 
with linear and constant terms included over those 
obtained by average loss formula excluding the linear 
and constant term, are marginal. 
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Find K for which D|^ is just greater than BAL 



C*=Con*Cln*X»C2n*x2-fF(N.,.K) 

:: * 


IF Cy <BEST >BEST =0^ i ST0 = x 



Fj^. =-1 > j = i + 1 to ISTEPS 



Compute individual generations Gjj 

i 

Confute transmission Losses for each Dj using 8^^^ 

i 

1 Compute load supplied for Dj 


I 


Find generation schedule of hydro plant using load 
levels and load supplied arrays and individual generations 

1 

r 

Compute 

Discharge 
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MEANINGS OF NOTATIONS AS USED IN THE PLOW CHART OIHER 
than those AIiREAD-J^ USED IN THE TEXT 


jp « Total nuiriber of plants 

qH * Nuinber of hydel plants 

:MDS « Minimum discharge specified 

jB = Upper bound of n^ plant 

^ = Lower bound of n^ plant 

n 

#Co *= Cost curve coefficients of n plant 
"Ol* In' 2n 


joad, 


3EST 


Bpsillon 


w Loss coefficient B-matrix 

= load level 

= i"^ demand in the derrand array 

= Cost of generation of plant at i^ demand 
level 

s Best solution so far obtained 
ss minimum lower bound 
= transmission loss at K load level 
» very small number 

= lowest value, of r (water value) say 0.0 at start 
» highest value of 7 say 1.0 at start. 
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